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Abstract. Proposed space-based gravitational-wave (GW) detectors such as DECIGO and 
BBO will detect ~ 10® neutron-star (NS) binaries and determine the luminosity distances to the 
binaries with high precision. Combining the luminosity distances with cosmologically-induced 
phase corrections on the GWs, cosmological expansion out to high redshift can be measured 
without the redshift determinations of host galaxies by electromagnetic observation and can be 
a unique probe for dark energy. This article is based on the results obtained in [l] where we 
investigated constraining power of the GW standard siren without redshift information on the 
equation of state of dark energy with future space-based GW detectors. We also compare the 
results with those obtained with other instruments and methods. 



1. Introduction 

Future space-based gravitational- wave (GW) detectors such as DECI- hertz Interferometer 
Gravitational-wave Observatory (DECIGO) [21[3] and Big-Bang Observer (BBO) [1] (see also [S] 
for updated information) are the most sensitive to GWs in 0.1 — IHz band and will aim at 
detecting the primordial GW background, the mergers of intermediate- mass black holes (BH), 
and a large number of neutron-star (NS) binaries in an inspiraling phase. 

It is known that the continuous GW signal from a compact-binary object provides a unique 
way to measure the luminosity distance to the source with high precision. Such binary sources 
are often referred to as the standard siren. With the redshift information determined by an 
electromagnetic (EM) follow-up observation, the standard siren can be an accurate tracer of 
the cosmic expansion [6]. The potential power of this method as a dark-energy probe has been 
investigated with ground- and space-based detector configurations for GWs [HI51I7HT7]. The 
most of the preceding works assume that the redshifts of all GW sources are known by EM 
observations. However, this assumption is rather strong and is too optimistic to be justified 
because the spectroscopic follow-up observation of galaxies is, in general, time-consuming, 
particularly at high redshifts [18]. In addition, not all galaxies can be observed due to intrinsic 
faintness, the absence of spectral features, limited sky coverage, and a limited redshift range 
(redshift desert). As a result, in a practical follow-up observation, the fraction of the binary 



sources whose redshifts are spectroscopically obtained is significantly reduced. From a rough 
estimate based on the number density of galaxies potentially observable and the number of 
galaxies to be observed in the future galaxy redshift surveys such as JDEM/WFIRST [19] and 
Euclid |20j, it turns out that the fraction of redshift identification is ~ 10~^ with uncertainty of 
about one order of magnitude [1]. This means that the GW sources with redshift information 
are quite rare unless we perform a large-scale follow-up campaign dedicated for GW events. 

Another approach is to measure the cosmic-expansion rate from GW observations alone 
without electromagnetically-estimated redshifts. As suggested in [2], the cosmic expansion 
affects not only the amplitude (luminosity distance) of GWs but also the phase, and it is possible 
to directly obtain information about the cosmic acceleration by accurately measuring the GW 
phase shift of a binary source at a certain redshift. Although the redshifts of the binaries are 
assumed to be determined by the EM follow-up observations in the previous works if 
we combine the independent information from the luminosity distance and the cosmological 
phase shift, we can measure the expansion history of the universe without any reference to 
the EM counterpart or host-galaxy identification. Furthermore, this method based on a purely 
GW observation enables us to compare the observational data with those obtained in other 
EM observations. Since proposed space-based GW detectors such as DECIGO and BBO would 
detect ~ 10^ NS binaries, they provide a novel opportunity to measure the property of dark 
energy without using the cosmic ladders. 

Our primary interest here is in the potential of the standard siren without any help of redshift 
information of the binary sources and how the sensitivity is improved by partially adding the 
redshift information of a fraction of NS binaries. We show the constraining power of the GW 
standard siren with/without redshift information on the equation of state (EOS) of dark energy. 
As GW detectors, we use future space-based detector, DECIGO and BBO, allowing their noise 
curves to be scaled appropriately and including confusion noises produced by astrophysical 
sources. Throughout the paper, we adopt units c = G = 1. 

2. Detector noise and astrophysical foregrounds 

In this section, we summarize the detector-noise curves of DECIGO and BBO and confusion 
noises from astrophysical GW foregrounds. 

DECIGO and BBO are the most sensitive to GWs in 0.1 - IHz band [2,4J. Using the 
currently-proposed design parameters 0, we obtained fitting formula for the sky-averaged noise 
curve of DECIGO single interferometer: 



where fc = 7.69 Hz. BBO has a slightly different noise shape from DECIGO, because of the 
different interferometer type and optical parameters. In the presence of WD confusion noise 
below ~ 0.1 Hz, the BBO sensitivity to a NS binary is nearly three times better in amplitude 
than that of DECIGO. 

DECIGO and BBO are sensitive to a large number of astrophysical GW sources, in particular 
NS, BH, and WD binaries. It is expected that the low-frequency side of the noise curve would 
be dominated by the GWs from extra-galactic population of WD binaries, which remains 
after subtracting individually identified signals and are stochastic in nature. According to 
the estimation by [22], the fitting formula for the power spectrum of WD confusion noise |23j 
(residual contribution after the subtraction) is given by 
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Figure 1. Noise curves of default and scaled DECIGO: r„ = 1 (blue), 1/2 (green), 1/3 (orange), 
and 1/5 (red) from the top to the bottom. These noise curves include the confusion noise from a 

number of WD binaries. The WD-binary foreground Sj^^ is shown with solid, black curve on 

the left side. The three diagonal dashed lines represent the NS-binary foreground \J~^^ before 
subtraction (or T^ns = 1) with no = 10~^, 10~^, 10"''' Mpc~^ yr~^ from the top to the bottom. 



There is also contribution from galactic WD binaries [23]. However, we verified that the galactic 
contribution is negligible above the frequency 5 x 10"^ Hz. 

Another astrophysical source that we have to take into account is the NS-binary foreground. 
According to [24J, the energy density of GWs from NS binaries per logarithmic frequency bin 
normalized by the critical energy density of the universe at present is written as 

Hq is the Hubble constant and is the chirp mass defined as Mc = rf^^Mt, together with the 
total mass Mt = nii + m2 and the symmetric mass ratio r] = mim2/Mf . h{z) is the NS merger 
rate per unit comoving volume per unit proper time at a redshift z. We adopt the following 
fitting form of the NS-NS merger rate given in h{z) = hQs{z), s{z) = 1 + 2z for z < 1, 
|(5 — z) for 1 < 2; < 5, for 5 < z. This function s{z) is estimated based on the star formation 
history inferred from the UV luminosity [2B]. The quantity tlq represents the merger rate at 
present. By assuming the flat ACDM universe with = 0.3 = 0.7) and the Hubble 
constant Hq = hj2 x 72kms~^ Mpc~^, we obtain the NS confusion-noise power spectrum 

5,f (/) ^ 1.55 . 10-.- [J^f' ( ,„_.„;:-3^-. ) . 

where we used the relation Sf^{f) = {3H^ /At:^ f)n^^{f) [27]. 

Including astrophysical contributions, a total-noise power spectrum is given by 

s,if) = sr\f) + 5WD(/) + snm^s . (3) 

Since our purpose here is to assess cosmology achieved by space-based detectors such as DECIGO 
and BBC and to make clear the requirement for the experimental design, we consider only 
DECIGO, but allowing the noise curve in Eq. ([1]) to vary by an overall factor r„ = 1, 1/2, 1/3, 
1/5 in amplitude. Namely, S'™'^*(/) = x «S'™q(/). Note that BBO noise power spectrum 



approximately corresponds to that of DECIGO with r„ = 1/3. The factor T^ns denotes a 
suppression factor due to the subtraction of individuahy identified NS binaries, whose value 
strongly depends on the detector noise curve. According to in which the authors have fitted 
the numerical result obtained by Yagi and Seto [28j, the suppression factor is T^ns = except 
for r„ = 1 case. When r„ = 1, some amount of the foreground residuals is left and gives 
T^NS ^ 4.6 X 10"'^. The contribution of each term in Eq. ([3]) is shown in Fig. [TJ 

3. Standard siren as a probe for dark energy 

3.1. Standard siren 

The continuous GW signal from a compact-binary object is often referred to as standard 
sirens [B]. For a single binary system, the Fourier transform of the GW waveform is expressed 
as a function of frequency / [271ES], 

Kf) = ;^M|/V-^/^e^*(^) , (4) 
dL\z) 

where is the luminosity distance, and the quantity = (1 + z)Mc is the redshifted chirp 
mass. Mc is the proper chirp mass defined in the source rest frame. The constant A is given by 
A = (\/6 TT'^/'^)^^, which includes the factor y^4/5 for a geometrical average over the inclination 
angle of a binary [30j. The function ^{f) represents the frequency-dependent phase arising from 
the orbital evolution, and at the order of the restricted 1.5 post-Newtonian (PN) approximation, 
it is given by p71[29] 



1 + y + ^^v) V-'^H^Mjf/' - IGTrrj-'/H^MJ) - ^X(.)M.(vr/M.)-«/3 

(5) 

where tc and (j)c are the time and phase at coalescence, respectively. The last term in the 
bracket is phase correction due to cosmic expansion (—4 PN order) [21l21j. where X{z) is defined 
as X{z) = ^[Hq — H{z)/{\ + z)] or equivalently expressed as [a(0) — a{z)]/2. Thus, X{z) > 
{X{z) < 0) corresponds to the accelerating (decelerating) universe. 

There are various methods of the standard siren. The simplest minimal one is to measure 
the amplitude of GWs, or the luminosity distance as a function of redshifts. However, seen from 
Eq. dH), GW observation itself cannot provide a source redshift, since the mass parameter we can 
determine is = (1 -|- z)Mc and the source redshift degenerates with the proper source mass 
Mc- So the redshift has to be determined from EM observation of host galaxies. This method 
as a tool to measure the cosmic expansion has been investigated by many authors [5t[7^ll4|. 
assuming that the source redshifts are known by spectroscopic follow-up observations. Even in 
the situation that the redshift of a source is not uniquely determined due to multiple candidates 
of a host galaxy in a detector error cube, one can infer the most likely cosmological model by 
statistically analyzing large samples of observed events [6l[31]. However, to do that, one needs a 
complete catalog of galaxy redshifts. This assumption is rather strong and is too optimistic to 
be justified, because of the smallness of the fraction of redshift identification (~ 10~^) pLj. 

To break the redshift degeneracy and measure the cosmic expansion by GW observation 
alone, recently some elaborate methods have been proposed. The redshift degeneracy could be 
broken if we know the EOS of a NS a priori [321 or if the observed chirp mass distribution of NS 
is sufficiently narrow |33j . We will discuss the sensitivities of several methods later in Sec. HI 

In this article, we investigate another method breaking the redshift degeneracy by the phase 
modulation of GWs due to the cosmic expansion, which is the last term including X{z) in 



Eq. ([5]) OEI]. Although in the preceding works the redshifts of the binaries are assumed to be 
determined by the EM follow-ups [SlEI]) the redshift information in principle is not mandatory 
if independent information from both the luminosity distance and the phase shift is available. In 
fact, as indicated by our results shown later, we can measure a cosmic-expansion rate only with 
GW observations without electromagnetically-estimated redshifts. This is a great advantage of 
this method, because a pure GW observation enables us to compare its observational data with 
those obtained in other EM observations. 

In the subsequent sections, we investigate the power of utilizing the phase correction without 
assuming any electromagnetically-estimated redshift, and give the figure of merits. Then we 
compare this method with the previous one, in which source redshifts are given, and combine 
both methods to maximize the accuracy of the estimation. 



3.2. Numerical estimation of measurement accuracy 

We consider a spatially flat universe with dark energy whose EOS is parameterized as 
w{z) = wq + Waz/{1 + z). The parameters to be determined from GW observations are wq, 
Wa, and Hq. For the analysis of error estimation, we adopt a fiducial set of cosmological 
parameters: wq = —1, Wa = 0, ftm = 0.3, Hq = 72kms~^ Mpc^^. 

In what follows, we calculate the estimation errors of binary parameters Oa- Mz, ij, tc, (pc, 
di, and X in the waveform of Eq. ^ with the Fisher information matrix. It is well known that 
the Fisher matrix analysis underestimates the parameter errors if SNR is not sufficiently large. 
Even default DECIGO (8 identical interferometers) is able to observe a NS binary at z = 1, 3, 
and 5 with SNR about 24, 31, and 67, respectively. Therefore, it is justified to use the Fisher 
matrix for the parameter estimation. The Fisher matrix for a single binary is given by 



where da denotes a derivative with respect to a parameter 6a. Since the 8 interferometers of 
DECIGO and BBO are assumed to be identical, the summation is reduced to just multiplying a 
factor 8. The noise power spectrum and GW signal are given by Eqs. ([3]) and dH), respectively. 
The lower cutoff frequency /mm is given by, /mm = 0.233 (IMq/Mz)^^^ (1 yr/T^obs)^^^ Hz, and 
the upper cutoff of the frequency is set to /max = 100 Hz which naturally arises from the upper 
cutoff of the noise curve. Given the numerically evaluated Fisher matrix, the marginalized l-cr 
error of a parameter, A^^! is estimated from the inverse Fisher matrix, — \/ ^"^00"' 

Once the estimation error of X is obtained, it is straightforward to calculate the estimation 
errors of the cosmological parameters. Since the error of dx, is much smaller than that of X 
(see Appendix of p]), we can replace observed di with the corresponding redshift in the fiducial 
cosmological model when we derive the measurement accuracies of cosmological parameters from 
X. Thus, for the simplicity of the analysis, we use X{z) instead of X[dL). Furthermore, we 
assume that the Hubble constant Hq is known a priori, and fix it to Hq = 72kms^^ Mpc^^, 
because the Hubble constant has been determined at a-few-percent level from the observation 
of nearby Cepheids and supernovae [35]. Thus, the free parameters of the Fisher matrix are wq, 
Wa, and Om- The Fisher matrix is given by 

_ r daX{z)d,X{z) dN{z) ^ dNjz) _ ^ Anr^jz) hjz) 
I al{z) dz dz "^""^^ H{z) 1 + z' 

where dN{z)/dz is the number of NS binaries in the redshift interval [z, z + dz] observed during 
Tobs [25]. r{z) is the comoving radial distance defined as r{z) = dL{z)/{l + z) and n is given in 
Sec. [2j Since the normalization of h is still uncertain, we adopt the most recent estimate [36], 
no = 10~^ Mpc^^ yr"^- Since the observation time Tobs is a crucial parameter, we suppose 3-yr, 
5-yr, and 10-yr observation. 








Figure 2. wq - Wa error ellipses marginalized over Vtm- The merger rate is fixed to 
no = 10~^ Mpc~^ yr~^. From the left to the right panels, the observation time is Syr, Syr, 
and lOyr. In each panel, the larger to smaller ellipses denote those with r„ = 1, 1/3, and 1/5, 
respectively. The dot at the center is our fiducial value: = —1 and Wa = 0. 



3.3. Results 

The measurement accuracies of wq and Wa marginalized over the other remaining parameters 
are shown in Fig. [5] for the fixed no varying the values of r„ and Tobs- The size of error 
ellipses change significantly depending on the parameters Tobs ^-nd r„ . For example, if we choose 
Tobs = 5yr, the errors Awo and l^Wa are 0.330 and 3.254 for r„ = 1, and 0.113 and 1.231 for 
fn = 1/3. For 10-yr observation, Aiuq and Aw^ reach at the levels of several xlO~^ and a few 
xlO~^, respectively. The more comprehensive parameter surveys are found in [IJ. 

To show the results more clearly, let us define the figure of merit (FoM) for dark energy, 
FoM = -y/detT^, where ^ab is the inverse matrix of (r^-'^)^^ with a,b = wq, Wa- This FoM 
is inversely proportional to the area of an error ellipse on the wq - Wa plane and ranges from 
~ 10"^ to ~ 10^. Since DECIGO and BBO will be launched in the late 2020s, the FoM should 
be compared with other future projects of an EM observation probing dark energy at that time: 
type-la supernovae, baryon acoustic oscillation, or weak lensing surveys. So typical criteria in 
FoM would be from 10 to 100, which correspond to future projects of stage III and IV in the 
dark energy task force [37j, respectively. To achieve these criteria, observation time longer than 

Q1 /Q 

5yr is preferable. Since the FoM rapidly improves being roughly proportional to T^^^ [T], the 
observation time is a crucial factor. Given 10-yr observation and the typical rate of binary 
mergers, the FoM is ~ 100 and is comparable to a stage-IV project in the dark energy task 
force. It should be emphasized that this method requires no redshift information of the GW 
sources and is completely independent of any EM observation. 



3.4-. Adding redshift information of binaries 

In the previous subsections, we did not assume any redshift information to measure the 
dark-energy parameters. Although it would be too ideal to assume that all binary redshifts are 
well determined by EM follow-up observations, some fraction of them will be determined by the 
follow-up observations. Then, we can optimize the constraints on the dark-energy parameters by 
using thus-determined redshifts. In this subsection, we combine z - di and di - X information 
and evaluate what fraction is needed to achieve FoM = 100, which is a typical value expected 
in future dark energy surveys [37] . 

The relevant Fisher matrix is defined with the redshift-identification fraction a by 

Jo ^dA^) dz Jq a\{z) dz 
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Figure 3. FoM as a function of a defined in Eq. ^ or A^^. The observation time and the merger 
rate are fixed to Tobs = 5yr and no = 10~® Mpc~'^ yr~^. The horizontal dashed line represents 
FoM = 100. The solid and dotted curves are the case with r„ = 1 and 1/3, respectively. 



where 0"^^ is an error in the luminosity distance. Strictly speaking, the variables di and X 
cannot be treated independently as in Eq. ([8|). However, the correlation coefficient between 
and X is of the order of 10~®. Hence, we can safely use the expression in Eq. ([8]). The fraction 
factor a would be a function of a redshift in a real galaxy survey, but we take it as a constant 
for simplicity of the analysis. 

We estimate the error size of the luminosity distance via the same procedure as in 
Sec. 13.21 but including possible systematic errors: weak-lensing magnification due to the 
matter inhomogeneities along the line of sight [38] and the Doppler effect due to the random 
peculiar velocities of binary sources [39]. These systematic errors to the luminosity distance are 
summarized as 



ri_ (1 + ^)^0.251 1-8 



(9) 
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where cJinst, ciens> and cjpv are induced by the instrumental noise, the lensing magnification, and 
the peculiar velocity of binaries, respectively. o"v,gai is the one-dimensional velocity dispersion of 
the galaxy and is set to be Cv^gai = 300km s~^. For the plot of the redshift dependence of each 
error term and the details, see [T6]. The important thing is that the lensing error dominates at 
almost all redshift range. 

Then using Eq. ([8]), we estimate the errors of wq and Wa and corresponding FoM as a 
function of a. We also calculate the corresponding number of NS binaries whose redshifts 
are identified. The result is shown in Fig. [3] for typical values of parameters, Tobs = 5yr and 
no = 10~^ Mpc~^ yr~^. At low a, of course, the z - di information plays no significant role 
in improving the FoM. But at a ~ a few xlO~®, z - di information begins to contribute to 
FoM, and enables the FoM to reach 100 at around a w a few xlO~^. As anticipated, z - di 
measurement is much more powerful than the di - X measurement because FoM|q,=o is much 
smaller than FoM|q,=i. Therefore, identifying the source redshifts strongly assists the di - X 
measurement. 

One may concern how much fraction a is feasible in the future galaxy redshift surveys. 
From our rough estimate [Ij based on the number density of galaxies potentially observable 



and the number of galaxies to be observed in the future galaxy redshift surveys such as 
JDEM/WFIRST [19] and EucHd [20], it turns out that the fraction of redshift identification 
is ~ 10~^ with uncertainty of about one order of magnitude. This means that the GW sources 
with redshifts are quite rare unless we perform a large-scale follow-up campaign dedicated for 
GW events. So it turns out that di - X measurement may be more sensitive than z - di 
measurement and that the standard sirens without redshifts guarantee minimally achievable 
FoM. 

4. Comparison with previous results with other detectors and methods 

Preceding works for other detectors are listed in Table [H in which only the detectors that 
can probe at cosmological distance are shown. As for DEICGO/BBO, our results should be 
compared with the work by Cutler and Holz [5]. They assume that one can identify a single 
host galaxy due to the good angular resolution of DECIGO/BBO and that the redshift of the 
host galaxy is always obtained by an EM follow-up observation. However, as estimated in [1], the 
assumption is too optimistic to be justified, because the redshift identification of host galaxies 
is time-consuming and inefficient and whose fraction would be at worst ~ 10~^ even for the 
future galaxy surveys such as JDEM/WFIRST and Euclid. Thus, our results are conservative 
because we do not require any redshift identification and any prior from other observations. We 
can conclude that the standard siren without source redshifts guarantee minimally achievable 
FoM, which can be improved by increasing the fraction of the galaxy redshift identification and 
combining other EM observations. 

One may also compare the results in this paper with those obtained from Einstein Telescope 
(ET). The total number of NS binaries observed by ET is comparable to DECIGO. However, 
ET has to rely on short gamma-ray bursts (GRB) to determine the source redshift so that the 
number of sources with redshifts is limited by short GRB directed toward us and is reduced to 
~ 10^ or smaller. Consequently, the sensitivity of DECIGO/BBO without redshift information 
is twice better than that of ET by Sathyaprakash et al. [12j. If the Planck prior of CMB is added 
as done by Zhao et al. [Hj, the size of the error is significantly reduced. On the other hand, 
the method proposed by Messenger and Read [32] does not assume source redshifts. However, 
to break the parameter degeneracy between the chirp mass and redshift of a NS binary, one 
needs to assume the EOS of a NS to be known. Furthermore, a tiny phase signal due to tidal 
deformation of NS (5 PN order) has to be detected. ET would be able to detect this signal, but 
the measurement accuracy of the redshift is not so good, Az/z ~ 0.1, and that of the luminosity 
distance also would be of the same order. To clarify its feasibility as a cosmological probe, 
further studies are needed. 

LISA has marginal angular resolution so as to identify a single host galaxy. Van Den Broeck 
et al. [13] considered the observation of a SMBH binary that fortunately accompanies a single 
host galaxy in its sky position error area. If only wq is a free parameter, the sensitivity is ~ 0.2. 
Although multiple sources improve this sensitivity, such fortunate events would be rare. To 
avoid this difficulty, a statistical analysis method has been developed [15]. This method utilizes 
all redshifts, which obtained a priori by EM follow-up observations, of possible host galaxies 
in a LISA error cube and determines the best fitting cosmological model for all GW sources. 
Petiteau et al. |15j performed the analysis with LISA, including the CMB prior, and obtained 
the accuracy of a few percents for wq, though other parameters are fixed in their analysis. The 
applications of this statistical method to ET and DECIGO/BBO are not straightforward. The 
problem is that the spectroscopic catalog of galaxies is far from a complete one for a large 
number of sources observed by these detectors. As mentioned above, the fraction of the redshift 
identification of host galaxies would be ~ 10~^. Therefore, it is interesting to ask which methods 
are more sensitive, performing the statistical analysis with an incomplete galaxy catalog or using 
all GW sources without redshift information as we have done in this paper. We should address 



Table 1. List of cosmological constraints obtained with other detectors and methods. Here we 
hst only the detectors that can probe at cosmological distance and works that take into account 
weak lensing error. SMBH and CMB refer to supermassive black hole and cosmic microwave 
background, respectively. 



reference 


detector 


(_rW 


number 


source redshift 


method to 






source 


of events 




determine z 


Van Den Broeck et al. [13j 


LISA 


SMBH 


1 


host galaxy 


identifying 












a single source 


Petiteau et al. [15] 


LISA 


SMBH 


~ 40 


host galaxy 


statistical 


Sathyaprakash et al. [T2] 


ET 


NS 


10^ 


EM counterpart 


identifying 










(short GRB) 


a single source 


Zhao et al. [H] 


ET 


NS 


10^ 


EM counterpart 


identifying 










(short GRB) 


a single source 


Messenger & Read [32] 


ET 


NS 


N/A 


unnecessary 


tidal effect 












on GW phase 


Cutler & Holz [5] 


DECIGO 


NS 


~ 10^ 


host galaxy 


identifying 




or BBO 








a single source 


Nishizawa et al. [1] 


DECIGO 


NS 




unnecessary 


cosmological 


(this work) 


or BBO 








phase drift 



reference 


free 

parameters 


fixed 

parameters 


prior 


sensitivity 


Van Den Broeck et al. [13j 


Wo 




none 




0.2 






^k, Ho 








Petiteau et al. [T5] 


Wo 


Wa, ^m, 


CMB 


/\wo ~ 


0.036 






Vlk, Ho 


(WMAP) 






Sathyaprakash et al. [l2j 


Wo, ^rn, ^k, 


Wa 


none 


/\wo ~ 


0.23 




Ho 










Zhao et al. [H] 


Wo, Wa, ^m. 


none 


CMB 


/\wo ^ 


0.099 




ilfc, Ho 




(Planck) 


AWa ~ 


0.302 


Messenger & Read [32] 


N/A 


N/A 


NS EOS 


Az/z ~ 0.1 










(for an 


event) 


Cutler & Holz [3] 


Wo, Wa, ^m. 




CMB 


Awo ~ 


0.01 




Ho 




(Planck) 


AWa ~ 


0.1 


Nishizawa et al. [1] 


Wo Wa, 


^k. Ho 


none 


Awo ~ 


0.113 


(this work) 








AWa ~ 


1.23 



this issue in a future work. 
5. Conclusion 

Proposed space-based GW detectors, DECIGO/BBO, can be a unique probe for dark energy. 
Using millions of NS binaries detected during the observation, we have estimated the sensitivity 
to the EOS parameters of the dark energy with/without identifying the redshifts of host galaxies. 



As a result, we found that the detectors without the redshift information have constraining 
power competitive to the future EM observations. This is a great advantage since the GW 
detector alone can probe for the cosmological expansion and enables us to compare the data 
of purely GW observation with those in other EM observations. With the help of the redshift 
information, FoM 100 corresponding to the stage III or IV of the dark energy task force is 
easily achieved with at most a few hundreds of sources. The standard sirens without redshifts 
guarantee minimally achievable FoM and strongly supports the feasibility of the space-based 
GW detectors as a dark energy probe. 
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